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Abstract 



Loop quantum cosmological methods are extended to homogeneous models in 
' diagonalized form. It is shown that the diagonalization leads to a simplification of 

the volume operator such that its spectrum can be determined explicitly. This allows 
| the calculation of composite operators, most importantly the Hamiltonian constraint. 

As an application the dynamics of the Bianchi I model is studied and it is shown 
^ ■ that its loop quantization is free of singularities. 

o 

a" 1 Introduction 

In applications so far, quantum gravity has mostly been analyzed in mini- or midi-superspace 
models which are obtained by ignoring an infinite number of degrees of freedom of the full 
theory. In particular, as a first approximation of our universe, which is homogeneous 
and isotropic at large scales, one often uses a model with a single gravitational degree of 
freedom, the scale factor or radius a of an isotropic universe To study more realis- 
tic models, and to see whether results of the highly symmetric approximation are robust, 
one has to break some of the symmetries. As a first step, one can break isotropy while 
retaining homogeneity. This provides a large class of models much more general than the 
isotropic ones which still have a finite number of degrees of freedom such that field theoretic 
complications can be avoided. 

The main question for quantum cosmology is the fate of the classical singularity. In 
isotropic loop quantum cosmology j2] the evolution has been shown to be singularity-free 
[3] , and to lead to additional unexpected quantum effects [HE]. To check the robustness 
of those results one has to see whether they persist in less symmetric, e.g. homogeneous, 
models. At the same time, homogeneous models show many different effects in their clas- 
sical approach to the singularity which may obtain quantum modifications. Homogeneous 
models also provide a hint at the behavior of the full theory: according to the BKL-scenario 
H different points on a space-like slice decouple from each other when one approaches the 
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singularity which suggests that it is enough to understand homogeneous models in order to 
see the fate of the classical singularity. Of course, the BKL-scenario itself is not established 
and so far it is purely classical. Still, the behavior of homogeneous models can provide 
suggestions as to the general approach to a classical singularity before we are able to study 
this issue in the full theory. 

Besides the reduction of degrees of freedom, the main simplification of isotropic loop 
quantum cosmology as compared to full quantum geometry is the fact that its volume 
operator is simple and its spectrum is known explicitly [7j. This is essential for direct 
computations since the volume operator plays an important role in constructing more 
complicated operators such as the Hamiltonian constraint [Hj. If one can deal with the vol- 
ume operator easily, one can also obtain an explicit Hamiltonian constraint equation which 
then can at least be implemented numerically. This has been essential in recent physical 
applications of isotropic loop quantum cosmology [3 13 IH EU 13 EJ U2] • Homogeneous 
models still have a finite number of degrees of freedom, but their volume operator 
corresponds to that of the full theory restricted to a spin network vertex with six incoming 
edges |7. . It is impossible to compute the spectrum of this operator explicitly since it would 
require a diagonalization of arbitrarily large matrices. 

However, this volume operator is that of a general homogeneous model which has six 
kinematical gravitational degrees of freedom (all components of a symmetric 3 x 3-matrix 
which represents the metric in a given point). In classical and standard quantum cosmolog- 
ical analyses one usually considers only diagonalized models with only the three diagonal 
components as kinematical degrees of freedom. Because the remaining components can be 
chosen to be zero by fixing gauge or symmetry transformations (depending on the model 
[T%]). diagonalized models are enough for physical applications. Diagonalizing a homo- 
geneous model is a restriction to a submodel since some degrees of freedom are ignored. 
In its physical interpretation this is different from a symmetry reduction, but the same 
techniques (developed in [T3]) apply. 

It turns out that this reduction suffices to simplify the volume operator resulting in an 
explicitly known spectrum. As in the isotropic case, this also requires to deal with gener- 
alized spin network states which do not appear as functions on a group manifold. One can 
represent those states in a convenient form such that an explicit analysis of homogeneous 
models is possible. The main application in the present paper is the result that the evolu- 
tion of homogeneous models is non-singular, which occurs in a way very similar to that in 
isotropic models. This is illustrated here in the Bianchi I model corresponding to a torus 
action of the symmetry group. 

Compared to isotropic models, the absence of a singularity in homogeneous models 
relies quite non-trivially on properties of quantum geometry. In isotropic models it has 
been seen to be important that quantum geometry is based on triad variables rather than 
the metric. This introduces a sign in the degree of freedom corresponding to the metric 
which is necessary to follow the evolution through the singularity. In homogeneous models 
we will observe the same effect, but it could also have been obtained if quantum geometry 
would be based on, e.g., the co-triad rather than the densitized triad. We will see that the 
removal of homogeneous singularities (contrary to isotropic singularities) would not happen 
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generically in a theory based on the co-triad. The reason is that the classical singularity 
lies in the interior of a minisuperspace built from the densitized triad, while it lies at the 
boundary (at infinity) of a minisuperspace based on the co-triad. 

After demonstrating how the diagonalization of a homogeneous model can be done 
in a connection formulation we will implement it at the quantum level. This requires a 
reduction of general homogeneous models along the lines of a symmetry reduction which 
is carried out explicitly. In this way we will obtain the kinematical framework whose 
structure is very similar to that of isotropic models, and derive the volume operator and 
its spectrum. This will allow us to compute the Hamiltonian constraint equation explicitly 
and to bring it into the form of a partial difference equation after transforming to the triad 
representation. As an application the classical singularity is shown to be absent. 

2 Diagonal Bianchi Class A Models 

Bianchi models are obtained by a symmetry reduction of general relativity with respect 
to a symmetry group S which acts freely and transitively on the space manifold E = S. 
Only class A models which have a symmetry group with structure constants fulfilling 
Cjj = are amenable to a Hamiltonian analysis. The action provides left-invariant 1-forms 
LO^ which yield the decomposition A l a = (f)\uj^ of an invariant connection with constant 
coefficients </>}. While the left-invariant 1-forms are regarded as background structure 
(provided by the symmetry group) with respect to which the reduction is being done, the 
components 0} describe the invariant connection in a reduced model. A corresponding 
decomposition of the densitized triad is given by Ef = pfXf where Xf are densitized left- 
invariant vector fields dual to the 1-forms. The symplectic structure of the reduced model 
is given by = ■yn5j5f where 7 is the Barbero-Immirzi parameter and k = 8nG the 

gravitational constant. (The role of a possible coordinate volume is discussed in |16j.) 

In this form, the invariant connection has nine independent components which, after 
subtracting three gauge degrees of freedom, yields six gauge-invariant components. This 
corresponds to the fact that a homogeneous model is described by <fr) in full generality 
without restricting to diagonal metrics. To perform the diagonalization we rewrite an 
invariant connection by using the polar decomposition <ft = Rip where R G 50(3) is a 
rotation matrix and ip a symmetric matrix. The symmetric matrix ip = Adiag(ci, C2, C3)A T 
can then be diagonalized introducing another rotation matrix A G 50(3). Defining O := 
RA G 5*0(3) and redefining ui'^ := uJ^Of we obtain 



Now we can clearly see that only three gauge-invariant degrees of freedom cj remain, which 
is just what we need for the diagonalization. (Note that we absorbed the rotation matrix 
O into the left-invariant 1-forms, which can be done in general without amounting to 
a reduction of the theory. In a reduced model, however, we have to fix the 1-forms in 
the beginning as a background structure, and then different forms differing by a rotation 
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would lead to different formulations of a model. This explains that the diagonalization 
leads to a reduction of the model, while any connection can be written in the form (JIJ 
without restriction if the uj are regarded as being free. Note also that the diffeomorphism 
constraint T>i = cj K (j)jpf c = cf K C(j)p( K ^5f vanishes identically in the diagonalized class A 
models whereas it generates inner automorphisms of the symmetry group in non-diagonal 
models.) 

In fact, the corresponding decomposition 

E a = p (K) A K x ,a (2) 

of the triad shows that the metric is diagonal. We have to allow negative values for the 
coefficients p 1 since the triad can have two different orientations while the orientation of 
A G SO (3) is fixed. The S'0(3)-matrix A contains purely gauge degrees of freedom and 
plays the same role as the matrix A in the isotropic decomposition [2] , even though it arose 
in a different way. The advantage of the diagonalization in the connection formulation 
is that gauge-invariant degrees of freedom in cj and gauge degrees of freedom in A are 
strictly separate (except for simple residual gauge transformations acting on the cj, and 
correspondingly on p 1 , which will be discussed in Sec. I3.1|) . As in the isotropic case, this 
will lead to a simplification of the volume operator since essentially SU (2)-calculations can 
be reduced to [/(l)-calculations. In the diagonal components the symplectic structure is 

{ C/ ,/} = 7^/- (3) 

Also the co-triad e l a = a\uj l a = a^k^u^ has the diagonal form with arbitrary real 
parameters aj G K. The relation between the triad and co-triad components follows from 
El = | dete|(e^) _1 and is given by 

P 1 = \a 2 a 3 \ sgn(ai) , p 2 = |a x a 3 | sgn(a 2 ) , p 3 = |a x a 2 | sgn(a 3 ) . (4) 

In these variables the volume is given by 



V = \j\\^ k ei JK p\p]p^\ = a/IpWI = laical • 

For the geometrical interpretation of the connection components we need the homoge- 
neous spin connection T l a = T^A)uj^. From the general formula 



we derive 



^L n J + ^L n K__^]_ n i} (5) 
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where the numbers n 1 specify the Bianchi model via Cj K = ejKLn IL = ^jx n • Only 
the Bianchi I model has vanishing spin connection; otherwise Tj depends on the triad 
components, contrary to isotropic models where spin connection components are constant. 
For diagonal models the Euclidean part of the Hamiltonian constraint (TBI E] reduces 

to 

= k- 1 det (^(eijkcfj^pi - ^jPH + titfjPjPk) 

= 2/iT 1 {{n l ci - c 2 c 3 )a! + (n 2 c 2 - CiC 3 )a 2 + (n 3 c 3 - CiC 2 )a 3 ) . (6) 

The full constraint H = -# (E) + P with 

P = -2(1 + 7 2 )*:- 1 det (e j ) " 1 K\j E- Ej , 

using the extrinsic curvature K] = 7 -1 (T} — 0}), is 

H = 2K- 1 (({c 2 r 3 + c 3 r 2 -r 2 r 3 )(i+ 1 - 2 )-n 1 Cl - 1 - 2 c 2 c 3 )a 1 

+ ((cir 3 + c 3 T 1 - rir 3 )(l + 7~ 2 ) - n 2 c 2 - ^ 2 cic 3 ) a 2 

+ ((cir 2 + c 2 T 1 - rir 2 ) (1 + 7 ~ 2 ) - n 3 c 3 - 7 ~ 2 cic 2 ) a 3 ) . (7) 

In terms of extrinsic curvature components Kj = \hi = 7 _1 (r/ — Cj), where the dot 
indicates a derivative with respect to proper time t corresponding to lapse AT = 1, it reads 

h = 2 K - 1 ((r 2 r 3 -n 1 r 1 )a 1 + (r 1 r 3 -n 2 r 2 )a 2 + (r^ - n 3 r 3 )a 3 

— \(a\a 2 a 3 + a 2 aia 3 + a 3 aia 2 )) (8) 

which yields the Hamiltonian constraint equation H+H raeAiter {ai 1 a 2 , a 3 ) = as a differential 
equation for the co-triad components aj. 

As special cases we obtain isotropic models if we set a\ = a 2 = a 3 =: a such that 

H + H mattcI (a) = -f^-^d 2 + 2I> + H mattei (a) = 

where T = for the flat model and r = | for the model with positive spatial curvature. 
In this case, the Hamiltonian constraint equation is the Friedmann equation. 

An example for an anisotropic model is the vacuum Bianchi I model where n 1 = and 
-^matter = 0. Therefore, the constraint equation reads 

dia 2 a 3 + a 2 aia 3 + a 3 aia 2 = . 

A simple solution is the one where the a/ are constant such that all terms in the constraint 
equation vanish separately; this solution is Minkowski space. All other solutions are given 
by the well-known Kasner solutions (201 a i ^ where the constant exponents ai have 
to satisfy Y2i a i = 1 = X]/ ?- These relations can be fulfilled only for —l<aj<l 
and, furthermore, one coefficient ai has to be negative while the other two are positive. 
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This leads to an evolution of the form of Kasner epochs where one direction of space is 
contracting while the other two are expanding. The total volume V = aia 2 a 3 oc t always 
increases monotonically. In these co-triad (or metric) variables aj the structure of the 
singularity at t = is quite complicated because one metric component diverges. This 
implies that the classical singularity lies at the boundary of minisuperspace. 

In the densitized triad variables p 1 on which the loop quantum cosmological descrip- 
tion is based the situation turns out to be different. Using the relation between co-triad 
and densitized triad components and the relations for the coefficients aj it follows that 
p 1 oc t 1- " 7 which has always a positive exponent. Thus, the triad components are always 
increasing and the singularity lies at p 1 = which is an interior point of minisuperspace. 
This will turn out to be crucial for the generic removal of the classical singularity. 

Other Bianchi models have a non-vanishing potential given by the spin connection 
components. This leads to possible changes in the expansion/contraction behavior of 
the co-triad components with intermediate regimes approximately described by Kasner 
epochs. Such a succession of different epochs can even be chaotic in certain models which 
renders the structure of the classical singularity even more complicated. In densitized 
triad variables, the situation looks again simpler since all components decrease during any 
Kasner epoch approaching the singularity. Only the component which decreases the fastest 
changes during the succession of different epochs. 

3 Quantum Theory of Diagonal Bianchi Class A Mod- 
els 

According to the general scheme of a symmetry reduction for diffeomorphism invariant 
quantum field theories of connections symmetric states are distributional states which 
are supported only on invariant connections with respect to the given symmetry group. 
For homogeneous models those states are in one-to-one correspondence with functions on 
the group manifold SU(2) 3 which can be visualized as spin network states associated with 
a graph containing a single vertex and three closed edges intersecting in the 6-vertex. For 
the diagonalization we have to perform an additional reduction to diagonal homogeneous 
connections ((TJ. Even though this is not a symmetry reduction, the same techniques apply. 
Consequently, diagonal homogeneous states are functions on the manifold 

{(exp(c 1 Air i ), exp^Ti), exp^A^)) :C(6R,Ag SO(3)} C SU(2) 3 (9) 

which is obtained by building (point) holonomies out of the connection components (this 
issue is discussed more carefully in [IS])- The diagonal homogeneous configuration space 
is by definition a submanifold of the general homogeneous configuration space SU(2) 3 , but 
as in the isotropic case [7j one can easily check that it is not a subgroup. This implies 
that the Peter-Weyl theorem cannot be used to find a generating set of functions, and 
states have to be expected not to be ordinary spin network states. In the isotropic case 
all possible functions of the lone gauge invariant connection component c were allowed as 
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gauge invariant states, not just even functions as would be expected by a naive reduction 
(corresponding to characters on a single copy of SU(2)). 



3.1 Diagonal Homogeneous Spin Network States and Basic Op- 



To derive all possible states for a given model one has to restrict the already known 
homogeneous states (or in general all states of the full theory) to the subspace of reduced 
connections, find an independent set, and complete it with respect to the measure on 
the space of reduced connections. The resulting set can be ort honor malized to derive 
a generalized spin network basis. In most cases this procedure can be simplified if a 
decomposition of reduced connections into pure gauge degrees of freedom and remaining 
components is known. Then, gauge invariant states can be found to be arbitrary functions 
of the remaining components subject possibly to residual gauge transformations which are 
usually discrete and easy to find. For instance, in the isotropic case an invariant connection 
has the form A l a = cA)u^ where the components A} are pure gauge. There are no residual 
gauge transformations on c (note that changing the sign of c is equivalent to changing the 
orientation of A, which is not a gauge transformation) such that gauge invariant states can 
be arbitrary functions of c without restrictions. This basic fact turns out to have important 
consequences even for the dynamics and physical applications: it implies that states are 
labeled by an integer quantum number n rather than a positive integer 2j + 1 as would be 
expected for an SU (2) -theory. In the end, this leads to the removal of classical singularities 
in isotropic models. 

We can now perform this analysis for diagonal homogeneous connections (0) which also 
allow a decomposition into pure gauge degrees of freedom A} and remaining components 
Cj. We start with noting the measure on the space of gauge invariant diagonal homoge- 
neous connections as a quotient of Since the relations between the three copies of 
SU(2) implied by the diagonalization only affect the gauge degrees of freedom A}, the 
gauge invariant measure is a direct product of three copies of the Haar measure on SU(2). 
Therefore, we use 



Gauge invariant states have to be functions of cj, but this time there are residual gauge 
transformations: it can be seen from Q that one can change the sign of two components 
Cj simultaneously while keeping the third one fixed (e.g., a gauge transformation with 
exp^AgTj) changes the sign of c\ and cz). Those are the only residual gauge transforma- 
tions such that gauge invariant states have to be invariant under changing two signs of the 

Cj. 

To write down a basis we use the isotropic states given by 



erators 



d/i(ci, c 2 , c 3 ) = (27r) 3 sin 2 (|ci) sin 2 (ic2)sin 2 (ic 3 )dcidc2dc3 . 



(10) 






7 



together with C~i( c ) — (v^2sin ~c) 1 and 



0(c) = cos^+|)c _ (i2) 



sm 



which are orthonormalized with respect to the measure d/i(c) = (27r) _1 sin 2 (ic)dc. Thus, all 
orthonormal states of the diagonal homogeneous model can be written as fj 1 (ci)gj 2 (c 2 )hj 3 (c 3 ) 
where /, g, h can be x or C an d ji is a non-negative half-integer for \ and a half-integer 
larger than —1 for (. However, only two different combinations are invariant under the 
residual gauge transformations, namely 

Xji,j2,j 3 ( C l^ C 2, c 3) : = Xji( c l)Xj 2 ( c 2)Xi 3 ( c 3) , (n,.j2,M C l> C 2, c 3) '■= Oi ( c l ) Cj 2 ( c 2 ) C73 ( C 3 ) 

(13) 

where the jj can be — \ only for £. As in the isotropic case, we obtain twice as many 
states as expected from positive spin labels. Such a doubling of the number of states has 
been seen to be necessary in the isotropic model because it yields eigenvalues for the triad 
operator corresponding to the two different triad orientations. To verify this in the diagonal 
homogeneous model we have to consider the quantization of the triad components p 1 . 

As in invariant vector field operators which quantize triad components in the full 
homogeneous model reduce to the self-adjoint derivative operators 

P / = -<7$(^ + 5«rt(ic/)) • (14) 

These operators, however, do not leave the space of states (jlHjl invariant corresponding to 
the fact that the triad components p 1 are not gauge invariant. As with connection compo- 
nents, we can change two signs of the p 1 simultaneously by a gauge transformation. Thus, 
only and sgnfj? 1 ^ 2 ^ 3 ) are gauge invariant, the latter quantity being the orientation. 
Still, for calculations it is convenient to work with the non-gauge invariant //-operators 
and to use their eigenstates which have to be non-gauge invariant. These eigenstates can 
easily be seen to be 

1^1,^2,^3} := \n>i) ® |n 2 ) ® |n 3 ) (15) 

where 

exp(|mc) 



(c\n) 



v^sinfic) 



■2 

are the isotropic p-eigenstates in the connection representation. Using Euler's formula for 
the exponential shows that the eigenstates (|15|) are not linear combinations of the gauge 
invariant states (fT3|) only. A gauge invariant state can be expanded in the states (fTK|) as 

\s) = ^2 s ni,n 2 ,n 3 \ni,n 2 ,n 3 ) (16) 

where the coefficients s niin2i „ 3 , which represent the state in the triad representation, have 
to fulfill 

(17) 
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The set of basic operators is completed by quantizing the connection components which 
has to be done via (point) holonomies hi = exp(cjA}Tj) = cos(|cj) + 2A}Tj sin(|c/) as 
multiplication operators. The matrix A is not gauge invariant, as a multiplication operator, 
and does not even leave the space of states invariant. Nevertheless, we do not need to 
extend our state space once more because, as shown in 9J, A drops out of gauge invariant 
operators even when holonomies are used in intermediate steps (as in commutators with 
the volume operator). We then need only the action of cos(|c/) and sin(|cj) which can 
readily be seen to be 

cos(±ci)|ni,n 2 ,n 3 ) = ±(|m + 1, n 2 , n 3 ) + |ni - 1, n 2 , n 3 )) (18) 
sin(|ci)|ni,n 2 ,n 3 ) = -§i(|ni + 1, n 2 , n 3 ) - |ni - 1, n 2 , n 3 )) (19) 

and correspondingly for c 2 and c 3 . 



3.2 Gauge Invariant Triad Operators 

As already discussed, the triad components are not completely gauge invariant, but only 
the functions \p ! \ and sgn(p 1 p 2 p 3 ) are. The corresponding operators must leave the gauge 
invariant states (fT5j) invariant. In fact, we have 

\P \Xh,h,h = rfpUl + ~i)Xhj2,h > \P \Cji,hJa = I^pUi + 2 )Cji,h,j3 

and 

fffxjuh^ = i (T^p) 3 (ii + |)(i2 + |)0'3 + |)0 1 ,i 2 j3 
fffCh,h,j3 = -*(7^) 3 (ii + |)0'2 + |)03+|)xiij2,j3- 
Simultaneous eigenstates of all gauge invariant triad operators are given by 

In terms of the states ()15|) the action reads 

\n 1 ,n 2 ,n 3 ) = ljlp\ni\ \nx,n 2 ,n 3 ) 

and 

sgn^Vp 3 )!^,^,^) = sgn(nin 2 n 3 )|ni,n 2 ,n 3 ) . 
Gauge invariant simultaneous eigenstates are 

n 2 , n 3 ) + \- n 1 , -n 2 , n 3 ) + | - n x , n 2 , -n 3 ) + K, -n 2 , -n 3 )) (20) 

= 2-3 (Cid^i-D.Ki^i-D.idnal-i) -isgn(n 1 n 2 n 3 )xi { | mhl)j i (|n2 |_ 1)j i (|n3 |_ 1) ) . 

From the triad operators we can immediately build the volume operator V = \/\p 1 p 2 p 3 \ 
(it can also be derived from the non-diagonal homogeneous volume operator [7j , which itself 
descends from the full one) with eigenvalues 



V\ni,n 2 ,n 3 ) = (^7/p) 2 y/\nin 2 n 3 \ \n 1 ,n 2 ,n 3 ) =: V ni)n2>ri3 \ni, n 2 , n 3 ) . (21) 
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To avoid confusion we note that the isotropic states \n) of [2| are not to be identified 
with states \n,n,n). Rather, according to the general scheme of a symmetry reduction 
an isotropic state would correspond to a distributional state in a homogeneous model and 
thus be non-normalizable. Therefore, the isotropic volume eigenvalues need not be given 

by K»,n,n- 

Using the basic operators we can construct more complicated ones like inverse metric 
operators along the lines of [HJ [SJ. Since there are no essential differences to the 
calculations in the isotropic model, we will not repeat details here. Those operators are 
needed when one quantizes the spin connection © or matter Hamiltonians. 

4 Dynamics of the Bianchi I Model 

So far we considered only kinematical aspects, which are independent of the particular type 
of Bianchi model. Specific properties only enter via the Hamiltonian constraint (JJj) (or the 
spin connection) which explicitly contains the structure constants of the symmetry group. 
Bianchi models usually have a non-vanishing spin connection which is necessary to yield 
non-zero intrinsic curvature. In the full theory such a statement would be meaningless 
because the spin connection does not have an invariant meaning. More precisely, it can be 
made arbitrarily small locally by choosing appropriate coordinates. In symmetric models, 
however, we are restricted to coordinate choices which respect the symmetry such that the 
spin connection does have a coordinate independent meaning in homogeneous models (or 
certain components in a less symmetric model). This also implies that it is impossible to 
ignore the effects of a possibly large spin connection whenever there is non-zero intrinsic 
curvature. As a consequence, even in classical regimes where, e.g., the extrinsic curva- 
ture is small, there can still be large intrinsic curvature contributions. Since this effect is 
a consequence of the symmetry and is not present in the full theory, it has to be dealt 
with appropriately. For instance, it affects the quantization of the Hamiltonian constraint 
(which has not been taken properly into account for the general quantization of the Hamil- 
tonian constraints in homogeneous models presented in (221) an d the interpretation of the 
semiclassical limit; this is discussed in more detail in [23] . 

The lone exception to the preceding statements is the Bianchi I model which has van- 
ishing intrinsic curvature and where the spin connection is always zero. Consequently, its 
quantization is simpler and also closer to that of the full theory [H] . In the present paper we 
restrict our attention to this model and in particular to the issue of the classical singularity 
displayed in the Kasner behavior. 

4.1 Hamiltonian Constraint 

From j22j we obtain the quantized Euclidean part of the Hamiltonian constraint in the 
form 

#( E ) = 4i( 1K 1*)- 1 ^e IJK tT(h I hjhJ% 1 hK[h] c \V]). 

UK 
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A further simplification of the Bianchi I model is that the Lorentzian constraint is related to 
the Euclidean part simply by H = j~ 2 H^ which we make use of here. Alternatively, one 
could derive the Lorentzian constraint as in the full theory jHj by expressing the extrinsic 
curvature as a commutator of the Euclidean part and the volume. Since this has been 
done in detail in the isotropic model without yielding essential differences to the simpler 
approach adopted here |2j, we will not perform the full analysis. In other Bianchi models 
the relation between the Lorentzian constraint and the Euclidean part is more complicated 
and one would have to quantize the additional terms (containing the spin connection) by 
hand, introducing more ambiguities |23j . 

Using diagonal (point) holonomies hi = cos(|cj) + 2A}Tj sin(|c/) in H := r )~ 2 H^\ we 
obtain 

H = 32ry~ 3 /ir 1 Zp 2 (sin(|ci) cos(|ci) sin(|c 2 ) cos(|c 2 ) 
x(sin(|c 3 )Fcos(±c 3 ) - cos(fc 3 )V sin(§c 3 )) 
+ sin(ic 3 ) cos(|c 3 ) sin(ici) cos(|ci)(sin(ic 2 )Ucos(ic2) — cos(|c 2 )V sin(~c 2 )) 
+ sin(|c 2 ) cos(ic 2 ) sin(ic 3 ) cos(|c 3 )(sin(|c 1 )V r cos(|ci) — cos(|ci)Vsin(|ci)) j 
with action 

H\ni, n 2 ,n 3 ) = 7~ 3 fiT 1 / P 2 ((V niiTl2)n3+1 - Ki 1)n2>n3 _i)(|ni + 2,n 2 + 2,n 3 ) 

-\m -2,n 2 + 2, n 3 ) - \n x + 2,n 2 - 2, n 3 ) + |ni - 2, n 2 - 2, n 3 )) 
+(V n 

1,712+1,713 l ,712 — 1,(13 

)(|n 1 + 2,n 2 ,n 3 + 2) 
-|ni - 2,n 2 ,n 3 + 2) - \m + 2,n 2 ,n 3 - 2) + |ni - 2,n 2 ,n 3 - 2)) 

1,7X2,113 v ni — 1,712,713 

)(|m,n 2 + 2,n 3 + 2) 
- |m, n 2 + 2, n 3 - 2) - |m, n 2 - 2, n 3 + 2) + |m, n 2 - 2, n 3 - 2))) . 

To derive an interpretation of the constraint equation (H + i^ ma tter)|s) = as an 
evolution equation, we transform to the triad representation (fTBjl where we have 

i.Hs^j ni ,712,713 K lp ((V^ ll _ 2n 2— 2 , 713+1 ^rti— 2,n2— 2,rt3 — l)^rti— 2,n2— 2,ns 



— (Kii+2,n 2 -2,n 3 - 


f 1 _ 




-2,n 2 


-2,713-1) 


"5 711 +2, 712— 2,713 


— (Kn -2, 712+2,713+1 — 


v ni . 


-2,712- 


+2,713-1) 


^711—2,712+2,713 


+ (Kii+2,7i 2 +2,n3- 


+1 _ 




-2,n 2 - 


+2,713-1) 


S 7ll +2,712+2,713 


+ (Kii-2,n 2 +l,(i 3 - 


-2 — 




-2,(12- 


-1,713-2) 


^711—2,712,713 — 2 


— (K ll +2,7i2+l,n 3 - 


-2 — 




"2,712 


-1,713-2) 


Srai+2,712,713— 2 


(^711-2,712 + 1,713- 


+2 — 


v ni . 


-2,712 


-1,713+2) 


^m— 2,712,713+2 


+ (Kii+2,n 2 +l,(i3- 


f 2 — 


Kll+2,(12- 


-1,713+2) 


^711+2,712,713+2 


+ (Kii+l,7i2-2,n 3 - 


-2 — 


Ki x - 


1,(12 


-2,713-2) 


"5(11,712—2,713—2 


_ (^711+1,712 + 2,713- 


-2 — 




-1,712+2,713-2) 


■8711,712+2,713— 2 


— (Kii+l,7i2-2,n 3 - 


f 2 — 


Kl- 


-1,712 


-2,713+2) 


■8)11,712 — 2,713+2 


+ (Kti+l,n 2 +2,7i3+2 - 




-1,712+2,713 + 2 


)Sni,n 2 +2,7i3+2 



(22) 
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Using the volume eigenvalues (jUJ), we obtain the evolution equation 



\n 3 + 2|a ^|n 2 + l|a - \n 2 - 1|H (Vi + 2| 2S ni+2 ,„ 2i n 3 +2 - 1^1 - 2|2S ni _ 2 ,n 2 ,n 3 +2) 

+ ( \n x + l|s - |ni - 1|5 ) ( |n 2 + 2| 2"s„ lin2+2jn3+2 - |n 2 - 2| 2"s ni ,n 2 -2,n 3 +2 



+ (jn 3 + l|a - |n 3 - l|aj Mn 2 - 2|a Mm - 2|2 Sni „ 2in2 _ 2in3 - |m + 2| a s ni+2in2 _ 2in3 
+ \n 2 + 2| 2 Mni + 2\ 2 S ni+ 2,n2+2,n a — \ n i — 2| 2 s m _ 2)n2 _|_ 2)n , 3 J J 

+|n 3 - 2|5 (f\n 2 + 1|* - |n 2 - 1|0 f|m - 2|5 Sril _ 2in2in3 _ 2 - |m + 2| 2"s ni+2in2 , n3 _ 2 ) 
+ (\m + l|s - - ^|n 2 - 2|^s Tll; „ 2 _ 2 ,n 3 -2 - |n 2 + 2 1 a s ni)n2+2in3 _ 2 jj 
= -2K(|7^) - 2^ inatter (n 1 ,n2,n3)s rilin2jn 3 (23) 

where ^matter can be any matter Hamiltonian which depends on the triad components 
ni and acts on the wave function s niin2jn3 (we suppressed the dependence of the wave 
function on matter fields). Equation (f2*Hj) is a partial difference equation for the wave 
function in the triad representation. For a dynamical interpretation we have to choose an 
internal time which in the case of the Bianchi I model can easily be done as one of the 
diagonal components of the triad, e.g., n.3. This is justified because classically p 3 is always 
monotonically expanding or contracting. With this internal time the constraint equation 
has the form of an evolution equation in a discrete time, as expected on general grounds 

EH- 

4.2 Absence of a Singularity 

The classical Kasner behavior suggests that the approach to the singularity of the Bianchi 
I model can well be described with the internal time p 3 (or any other diagonal triad 
component). In the loop quantization this leads to the discrete internal time n 3 which we 
already chose to interpret Equation (J23)) as a time evolution equation for the wave function 
s ni ,n 2 ,n 3 - The other labels n\ and n 2 can be regarded as independent degrees of freedom 
evolving in the time n 3 . While the time variable n 3 is unrestricted in the range from —00 
to +00, for independent values of the wave function we have to consider only non-negative 
ni, n 2 > 0. Remaining values of the wave function are then given by the relations 



-ni,n 2 ,n 3 ^rti,n 2 ,— n 3 ; °ni,— n 2 ,n 3 ^ni,n 2 ,-n 3 



(24) 



which follow from gauge invariance. In other words, s can be seen as a function on a 
discrete space {(ni, n 2 , n 3 ) : ni, n 2 e No, n 3 6 Z} subject to the identifications (ni, 0, n 3 ) = 
(rii, 0, — n 3 ) and (0, n 2 , n 3 ) = (0, n 2 , — n 3 ). Certainly, any one of the discrete coordinates can 
be chosen to take negative values; however, it is most convenient to choose that one which 
is being used as internal time. As in isotropic loop quantum cosmology, the difference to 
the classical and Wheeler-DeWitt approach is that there is a "doubling" of minisuperspace 
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which comes from the two possible orientations of a triad. This can be obtained in any 
triad formulation, but the non-trivial fact is whether or not the evolution can be extended 
meaningfully and uniquely through the boundary corresponding to the classical singularity. 

To study this question we use the evolution equation (J2HJ). We will interpret this 
equation as an evolution equation in the internal time 713, but it is important to note 
that the presence of a singularity is independent of the choice of internal time. The 
question is whether or not the wave function can be extended uniquely throughout all of 
minisuperspace even crossing singular hypersurfaces. For an intuitive picture one usually 
selects a function on minisuperspace, most often the volume, as a time coordinate and 
views solutions in a semiclassical regime as representing wave packets evolving in this 
internal time. But this is only necessary for an intuitive interpretation, not for a proof of 
absence of singularities. In particular, one can choose a time coordinate different from the 
conventional volume. Here we use the integer n 3 , corresponding to the triad component p 3 , 
because it leads to a difference equation, while the volume does not have equally spaced 
eigenvalues and thus does not give a simple evolution equation. Physical results, e.g. 
expectation values of observables, can always be translated between different internal time 
pictures. 

As a recurrence relation, equation (J2*3*)) does not only require initial conditions at slices of 
constant n 3 but also boundary data for small values of n\ and n%. We will first assume that 
the correct amount of boundary conditions is given such that we can solve the evolution 
equation for positive n 3 and negative n 3 , respectively. The issue of initial and boundary 
conditions will be discussed later. 

Starting with initial data at large positive times n 3 and evolving backwards, equation 
can be used as a recurrence relation as long as the coefficient \J\n% — 2| of the lowest 
order part, which contains the wave function with label n 3 — 2, is non-zero. This is the case 
only until n 3 = 2 such that our initial and boundary data uniquely fix the wave function for 
all n 3 > 0. The values at singular states, identified with n 3 = or n\ = or ri2 = 0, cannot 
be determined which looks like a breakdown of the evolution at the classical singularity. 
However, it turns out, in exactly the same way as in the isotropic case, that we can evolve 
through the singularity: the values of the wave function at the classical singularity remain 
undetermined, but since they completely drop out of the evolution equation they are not 
needed to find the wave function at negative times n 3 . The initial data chosen at positive 
n 3 , together with boundary data, uniquely determine the wave function at positive and 
negative n 3 and the evolution through the classical singularity is perfectly well-defined. 
Here we also used the fact that matter Hamiltonians quantized with quantum geometrical 
methods are zero at the classical singularity such that also the right hand side of 
(|2*3*|) does not contribute values of the wave function at the classical singularity. Again, 
it is essential that the orientation of the triad provides us with twice as many states as 
expected in a metric formulation. Thereby we naturally obtain negative values of our 
internal time. Still, just using variables which can take negative values is not enough 
since, most importantly, the evolution must remain well-defined. This is a non-trivial fact 
which in this case relies in part on the discreteness of the triad minisuperspace. As already 
alluded to earlier, another non-trivial fact (which is not present in the isotropic case) is 
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that the classical singularity lies in the interior of minisuperspace in the densitized triad 
variables which we have to use in quantum geometry, while it would be at the infinite 
boundary for, e.g., co-triad variables. It would be harder to see in general if a singularity 
at the boundary of minisuperspace can be eliminated with the present mechanism. 

We also note that choosing a different factor ordering of the constraint would lead to a 
singularity since some values of the wave function could not be determined from the initial 
data but now would not drop out completely and thus be necessary to determine other 
values. The ordering we have to use here for a non-singular evolution, corresponding to 
ordering the triads to the right, is the same as in the isotropic case. 

In order to discuss the boundary conditions we first transform the evolution equation 
by introducing 

and defining 

d(m):={ \l\ l + \ m ' l \-\l\ l -\ m ' l \ form^O (2g) 
[ for m = 

Since the volume eigenvalues V niin2tri3 vanish for n\ ri 2 ^ 3 = 0, the function t has to fulfill 

for all mi, 7712, J713. In addition, t also has to fulfill the gauge invariance conditions (|24j) 
such that we can restrict again to non-negative mi, m 2 . Equation (|23|) then turns into 



d(m 2 ){t 

mi+l,Tri2 ,7713+1 ^mi— l,fH2,wi3+l) d\TTli) (^mi,W2+l)TTi3+l ^mi ,m,2 — 1,?7J3+1 ) 

~\~d\Tn%) \tfn\— l,jna— l,m3 ^mi+l,m2— l)Ti3 "I - imi+l,m2+l,m3 tmi— 1,1712+1,1713) 

-\-d\Tfl2) (t JTtl _i jm2) }Tt3 — 1 ^mi+l,m2,m3-l) "I - d(wi]J (^mi,m2 — l,m3— 1 ^mi,m2+l,JTO3 — l) 

-7 3 K/p p matte r (mi, m 2 ,m 3 ) t 

mi ,ni2 ,m3 

(27) 



where 



Pmatter(mi, m 2 , m 3 ) = H mattei ( mi , m 2 , m 3 )/V mi>m2jm3 (28) 

for non-zero mi, m 2 , m 3 and zero otherwise, is the matter energy density operator. 

Starting with initial data at two slices of constant m 3 yields via recurrence the combi- 
nation 

d[Tn%) (^mi+l,m,2,m3+l tm\ — l,m 2 ,TO3+l) ~t" d\TTlx) (^mi,m2+l,JTi3+l ^mi,m2- l,W3+l) 

of values of the wave function at the next step m 3 + 1. For this to yield unique values 
V,m2,m 3 +i we need to specify boundary values at low values of mi and m 2 . For mi = 
we obtain the values of 

(f(m 2 ) (^l,m2,Jn3+l t— l,m2,m,3+l) <^(m 2 ) (^I,m2,m3+1 tl,rri2,— 1713 — l) 

from which we can find £i, TO2 , m3 at negative m 3 provided that the values £i, TO2 ,m 3 f° r positive 
m 3 are given by boundary conditions. The value mi = 1 in the recurrence relation yields 

<i(m 2 )t 2j . m2im3+ i 

,m 2 +l,m3+l ^l,vn,2 — l,m.3+l ) 
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(using to,m 2 ,m 3 +i = 0) which allows us to determine t 2im2>m , 3 +i uniquely. All other values of 
t for mi > 2 are fixed in the same way, and a similar analysis holds for varying m 2 . This 
shows that we need to specify the values 

ti,m 2 ,m 3 and £ mi ,i, m3 for all m 1 ,m 2 , m 3 > (29) 

as boundary data together with initial values at fixed 7773. 

As for the initial values, one would expect to need the wave function at two slices of 
constant 7773 because the evolution equation is of second order. However, by definition 
t mil m 2) m 3 has to vanish for 7773 = which is equivalent to the fact that s ni>n2) o drops out of 
the evolution equation. Thus, arbitrary initial values at fixed and 777.3 + 1 wm n °t be 
consistent with this condition and the evolution to 7773 = will yield linear relations for the 
initial data. Consequently, only the values at one slice of fixed 777.3 are f ree an d the number 
of initial conditions is only half the expected amount, as in the isotropic case. Unlike the 
isotropic situation, however, this does not suffice to specify the gravitational part of the 
wave function uniquely up to norm. Further conditions can be expected to arise from the 
requirement of a pre-classical wave function at large volume jl] where the wave function 
should not oscillate on Planck scales. This also restricts the original wave function s ni)n2jn3 
at odd values of rij, which are not determined by £ mi , m2i m 3 . Such an analysis, however, 
is more complicated for a partial difference equation and will be pursued elsewhere. Figs. 
^ and 121 show two different impressions of a solution to the Hamiltonian constraint of 
the Bianchi I LRS model which represents a wave packet moving through the classical 
singularity. Details of this model and its classical and Wheeler-DeWitt behavior, as well 
as more informations about the figures, can be found in the Appendix. 



4.3 Flat Space-Time 

It is usually hard to find exact solutions of a partial difference equation as (|2*Tj) with non- 
constant coefficients, even if it is linear. In our case it is easy to see that t mi , m2 , m3 = mi, 

tmi,m2,m:j 

= m 2 , and t mi , m2 , m3 = 7773 are exact solutions in the vacuum case. However, they 
do not satisfy the conditions (f2*U|) and thus are not allowed. 

An approximate solution at large volume can be seen to be of the form t mitTn2iTn3 = 
777i77727773. It does fulfill the conditions (J26|) . but this is not relevant since the solution is 
valid only approximately at large mi, m 2 , 7773. To see that this is an approximate so- 
lution we expand d{m) = \rri~ 1 + 0(m -2 ) such that the left hand side of (}2Tj) becomes 
the constant 6 up to terms of the order mj 1 . This corresponds to a negative matter en- 
ergy density p mat ter ~ — 67 _3 K _1 /p 2 V r ~ 2 which decreases like the inverse second power of 
the volume and thus would be negligible compared to any other standard matter compo- 
nent at large volume. In fact, the solution t mi)m2)m3 = mim 2 m 3 , which corresponds to 
s ni,n2,n 3 °t ^S n ( n i n 2 n 3) \/l n i n 2^3| oc sgn(77i?7 2 7i3)Ki lin2in3 , represents flat space in the sense 
that all terms in the constraint Cic 2 a3 + c 2 C3ai + 030103 = which follows from (J2J) vanish 
separately. At large volume a quantization of 03 = >y\p 1 p 2 /p 3 \ sgn(p 3 ) must have approxi- 
mate eigenvalues of the form a/ 1 77^/773! sgn(n3) such that in the triad representation we 
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Figure 1: Contour lines for a wave packet of the Bianchi I LRS model moving through 
the classical singularity (the diagonal line where the wave function drops to zero). The 
lower right boundary corresponds to the boundary m = of the minisuperspace, while the 
upper left boundary has been introduced for numerical purposes. (Not all values of the 
wave function are plotted; see the Appendix for further explanation.) 




Figure 2: Closer look at the classical singularity. Lattice points are values of the discrete 
wave function t mn . 
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have for our approximate solution (a 3 s) ni ^ n , 2t n 3 oc nin 2 . In the connection representation, 
this corresponds to the state 

a 3 |s> = E n i n 2 \ n i, n 2,n 3 ) oc S(ci)S(c 2 )0(c 3 ) 

ni,n 2 ,n 3 

where 5(c) := X^(2j+l)Xj( c ) denotes the ^-distribution on SU (2) and also 0(c) := J2j Cj( c ) 
is to be interpreted in a distributional sense. Therefore, one would expect that 63 |s) is 
annihilated by a quantization of C\c 2 as a multiplication operator, which appears in the 
constraint. However, as already seen, the annihilation is not complete which leads to the 
matter contribution decreasing as the inverse second power of the volume. To understand 
this fact we have to take into account that the configuration space (for a single diagonal 
component) differs from SU(2) precisely in the singularity structure of the point c = 
which is a fixed point of the gauge transformations. Correspondingly, the state space for a 
single diagonal component is not given simply by gauge invariant functions on SU(2), i.e., 
only the characters X31 but by all possible functions generated by Xj an d £r Since Q is 
not regular in c = for any j, an ordinary 5-distribution on this space of states cannot be 
well-defined. Instead, the distribution 5 introduced above appears which can be regarded 
as a regularization of the ill-defined naive (^-distribution. 

For regular test functions /(c), 5 in fact has the expected action, as can be checked in 
the usual way: 

(27T)- 1 J ~8(c)f(c) sin 2 (|c)dc = (27T)- 1 E I (2j + 1) sin ((j + \)c) sin(§c)/(c)dc 

= ^ £/ o cos((j + i)c)-(/(c)sin(ic))dc 



j Jo 

»2w 



l>ZTT 

= 2tt^ 1 V / cos(nc)(f (2c)sinc+ |/(2c) cosc)dc 
= 2(/'(2c) sine + \f(2c) cose) | c=0 

r-2-K 

-it- 1 / (f'(2c) sine + \f(2c) cos c)dc 
Jo 

= /(0) 

where we used the U(l) 5-distribution S(t) = (2tt)- 1 Y. n ei eM = ( 27r ) _1 ( 2 EneN cos(nf)+l) 
and periodicity of / and /'. If / is not regular in c = 0, however, i.e., of the form 
f(c) = g(c) sin(|c) _1 with a periodic function g which is regular and non-zero in c = 0, we 
have 

(27T)- 1 j ~8(c)f(c) sin 2 (ic)dc = (27T)- 1 J2 I (2J + 1) sin {(j + \)c) g(c)dc 

= 2tt- 1 J2 cos(nc)g'(c)dc=2g'(0). 



17 



Thus, the regularized ^-distribution is 



m = 4 (sin(|c)/(c)) u = { /jo) ;[ ; ( ^ ( • (*» 

Since the (^-distribution is regularized in order to take into account non-regular test 
functions, it is not annihilated when multiplying with sin(~c) or c. Instead, we have 

//*27T 
sin(i C )5(c)/(c)sin 2 (i C )dc = 27T 1 / cos(nc)(/' (2c) sin 2 c + /(2c) sine cos c)dc 

■A) „ cKJ 



nGN 

if / regular 



but 



(27r)- i y'sin(ic)5(c)/(c)sin 2 (ic)dc = {2-n)- 1 J 5(c)g(c) sin 2 (§c)dc 

= g(0) otherwise. 

In fact, it is easy to check that we obtain the distribution defined above when we multiply 
with sin(^c): 



0[/] = (27T)- 1 / 6(c)/(c)sin 2 (±c)dc = (2tt)- 1 ^ f cos(±nc)/(c) sin(±c)dc 
= (/(2c)sinc)| c=0 = 



if / regular 
S(O) if /(c) = 2(c) sin(ic)- 1 • 

Thus, sin(|c)5 = and we have to multiply a second time to get zero, sin(|c)0 = 0. 

In the isotropic model we have the same structure of the singularity c = in the 
configuration space, but there the constraint is quadratic in the single gauge invariant 
connection component c which annihilates the regularized ^-distribution. Therefore, there 
is no remaining matter contribution in the isotropic case. This demonstrates how certain 
features of a theory can be lost when additional symmetries are introduced, and also 
that working with gauge-invariant SC/(2)-states (fTTjl and (fT2"j) can lead to results different 
from those obtained with £7(l)-states e mc which would appear in a gauge-fixed theory. In 
this example, however, the physical relevance is not clear. Furthermore, it is an effect 
of the singularity structure in cj = which we can probe here since we are looking for 
solutions of flat space-time which are cj = classically. For other solutions the set cj = 0, 
which is of measure zero in the minisuperspace, is not expected to contribute significantly. 
However, such solutions corresponding to the classical Kasner behavior are harder to find 
even approximately, and one has to resort to a numerical analysis. Some difficulties to be 
faced are discussed in the Appendix. 
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5 Discussion 



In this paper we developed loop quantum cosmological techniques for homogeneous cos- 
mological models. The calculations have been simplified considerably by diagonalizing the 
canonical degrees of freedom, a process which retains all the relevant physical information. 
Most importantly, the volume operator is simpler and it is possible to find its spectrum 
explicitly. Having a simple volume operator then also leads to explicit expressions such 
as inverse metric operators and the Hamiltonian constraint operator which are necessary 
to study the classical singularity. Extending the results for isotropic models, we now have 
practical techniques for a larger class of systems given by all Bianchi class A models and 
their locally rotationally symmetric submodels (which have a single rotational axis). Those 
models are dynamically non-trivial even in the vacuum case such that they allow us to study 
quantum geometry effects independently of the matter coupling. 

Here we mainly applied the techniques to the Hamiltonian constraint of the Bianchi I 
model. It turned out that the classical singularity is absent in loop quantum cosmology 
in very much the same way as in the isotropic case. In particular, it is the same factor 
ordering of the constraint operator which is selected by the requirement of a non-singular 
evolution. New features compared to the isotropic case also emerged: it has been seen to 
be essential that quantum geometry is based on densitized triad variables, and the issue of 
initial conditions is more complex. One has to choose not just initial conditions but also 
boundary conditions on minisuperspace, and the issue of pre-classicality is more compli- 
cated to analyze. It may turn out that pre-classical wave functions (up to norm) fulfilling 
the consistency conditions are not unique, but still the freedom is reduced compared to 
the Wheeler-DeWitt approach. 

The analysis has also been simplified by choosing a diagonal component of the triad 
as internal time, rather than the volume. Even though the volume spectrum is known 
explicitly, its use as internal time would lead to a more complicated evolution equation 
since it is not equally spaced. For the Bianchi I model with its classical Kasner behav- 
ior, a diagonal component as internal time is justified, and it is well suited to study the 
approach to the classical singularity. In other Bianchi models the behavior is certainly 
more complicated: one expects a succession of different Kasner epochs. While in all those 
epochs the triad components p 1 are all decreasing toward the singularity, it is not always 
clear what happens during transitions between different epochs. If the behavior of the 
triad components remains monotonic during the transitions, one can use one of them as 
internal time and follow the evolution. Note that the result of an absence of a singularity 
will hold true in any case since it is a statement about the possibility of extending the wave 
function throughout all of minisuperspace uniquely; choosing an internal time is then only 
necessary to develop an evolution picture. 

Particularly interesting is the Bianchi IX model which has been suspected to behave 
chaotically. Here one can check if the modified approach to classical singularities in quan- 
tum geometry, which has been observed in the isotropic case in j3] and been seen to lead 
to inflation, implies a different behavior. In isotropic models, the origin of a modified 
approach comes from quantum geometry effects in the matter Hamiltonian while the grav- 
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itational part of the effective classical constraint equation has been left unmodified in a 
first approximation. In Bianchi models other than type I one expects a similar modifica- 
tion even in the vacuum case because the spin connection @ is non-zero and depends on 
triad components. This has to be dealt with properly when quantizing the Hamiltonian 
constraint and opens the door for modifications in effective classical equations of motion. 
While this does not modify the kinetic part (i.e., the terms containing time derivatives in 
(JHJ)) of the Hamiltonian constraint in a first approximation, it leads to a modified potential 
containing the spin connection at small volume. One can then expect that the series of 
transitions between Kasner epochs, which can be seen as reflections of the scale factors 
on the potential well, will change since the potential well breaks down if the volume is 
small enough. This may stop the reflections altogether and lead to a non-chaotic behavior. 
Then, a single Kasner epoch would be responsible for the transition through the classical 
singularity, which would not look much different from the Bianchi I behavior close to the 
singularity. These issues will be discussed in more detail elsewhere. 

The Bianchi IX model also provides further clues for the full theory. It has been 
speculated that the approach to inhomogeneous singularities is described by the BKL 
scenario according to which different points on a space-like slice decouple from each 
other close to a singularity (technically, only time derivatives remain in an asymptotic 
expansion of Einstein's equations close to a singularity). If this scenario is true and can be 
used in the quantum theory, it suggests that a non-singular Bianchi IX behavior implies 
non-singular behavior in the full theory. However, since even the classical validity of the 
BKL scenario is far from being established, such a conclusion would be premature at 
present. 

An alternative route is to study less symmetric models and see if what we have learned 
from more symmetric ones remains true. In this paper we presented the first step, going 
from isotropic models to anisotropic but still homogeneous ones. It turned out that essential 
results did not change, in particular the evolution remains non-singular. Even though the 
mechanism is similar to that in isotropic models, there are non-trivial aspects which cannot 
be seen in an isotropic model. It could well have happened that there is a singularity in 
homogeneous models but not in isotropic ones, e.g., if one used a theory based on the co- 
triad rather than the densitized triad as fundamental variable. That quantum geometry so 
far has provided precisely those properties one would need to have non-singular behavior 
is encouraging, but it also indicates that there may still be further surprises when we go 
one step further by studying inhomogeneous midisuperspace models. 
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Appendix 



A The Bianchi I LRS Model 

Probably the most simple anisotropic model is the Bianchi I LRS model without matter, 
which is defined as the subclass of Bianchi I geometries admitting one spatial rotational 
symmetry axis. Therefore, two of the diagonal components of the connection as well as of 
the triad, e.g. the first two for definitiveness, have to equal each other and only two degrees 
of freedom are left which we choose to be (c,p c ) and (A,pa) embedded into the Bianchi I 
model by 

(A, c) i-> (ci, c 2 , c 3 ) = (A, A, c) , (pa,Pc) ^ (p\p 2 ,P 3 ) = (pa,Pa,Pc) ■ 
The symplectic structure derived from the embedding follows from (JHJ): 

{A,p A } = ^K , { c ,p c } = 1 K. 
For the LRS model the constraint (J7|) reduces to 

H = -27~ 2 K~ x \p e \~% A(2cp c + Ap A ) oc Ap A + 2cp c 
which implies the equations of motion 

Pa = \pa , A = -\A , p c = 2p c , c = -2c 

solved by 

Pa(t) oc \fr , A(t) oc , p c (r) oc r 2 , c(r) oc r -2 . 

i 

Using p c as internal time we can write the evolution of p A as Pa{Pc) oc pi, which coincides 
with the Kasner behavior restricted to the LRS case: If two of the Kasner coefficients are 
equal, we have 2«i + «3 = 1 = 2a\ + a\ implying a% = | and a 3 = — |. In the Kasner 

1 4 

time coordinate t we then have p A oc ts, p c oc £3 according to (j4J), which also leads to 

1 

Pa(Pc) oc pi. 

The Wheeler-DeWitt equation in the triad representation for the Hamiltonian ApA + 
2cp c can be written in the form 

^Wa Va + 2 w) ^' Pc) = 

where a factor of p c has been absorbed in the wave function 4>{pa,Pc) '■= Pc^iPAiPc)- 
Fig. El shows a solution which represents a wave packet following the classical trajectory 

Pa{Pc) oc pi in minisuperspace, hitting the classical singularity. 

In the loop quantization we have only two integer labels m and n, where m represents 
the value of the two equal triad components pa and n the third one, p c . Residual gauge 
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Figure 3: A solution to the Wheeler-DeWitt equation of the Bianchi I LRS model. When 
the wave packet approaches the classical singularity (upper right corner), its amplitude 
diverges. The minisuperspace boundary = is at the right. 

transformations can now be seen to change the sign of m such that we can get unique 
representatives of gauge equivalence classes by considering only non-negative values of m. 
The label n, on the other hand, is allowed to take any integer value and thus is best suited 
to describe the transition through the classical singularity. We will use it as internal time, 
the sign of which, sgn(n), is again the orientation. 

For the sake of simplicity we write down directly the evolution equation for a wave 
function t m ^ n in the form (}2*Tj) which now becomes 

t n +l,m—l ) + d(n)(t 

^n— l,m+l ) 

= 2d2{m)t n+ i^ m + d(n)(t ntm+ i — £ n ,m-i) — 2d2(m)t n -i tm = . 

Here, we defined i ntTn := i n ,m+i — t n ,m-i an d dzim) := m~ l for m ^ and ^(0) := 
(it replaces d(m) of (|23j) since now two labels in the volume eigenvalues are equal to each 
other) . 

We now have the evolution equation in the form of a partial difference equation in two 
independent variables, which is straightforward to implement numerically. It is, however, 
a non-trivial problem to find numerical solutions which are close to pre-classical ones, i.e., 
which vary significantly only on large scales. We have many parameters for a solution, 
initial values at one n-slice and boundary values at small m. To find a solution which 
satisfies the consistency condition (to,m = 0), we can choose the initial values at the n-slice 
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and the boundary values freely. A complete initial value problem would be specified by 
also choosing values at the next n-slice since the evolution equation is of second order in 
n. We can find the values of the second slice by first leaving them as free parameters, 
solving the evolution equation up to n = where we obtain £ ,m as linear functions of the 
free initial values, and solving the set of linear equations to,m = 0. The result is a solution 
which by construction satisfies the consistency condition. In practice, we have to face the 
following problems when we want to find a numerical solution which is also pre-classical: 

• There is no guarantee that the computed consistent initial data at the second n-slice 
are close to the ones chosen at the first slice (which would be necessary for pre- 
classicality). One can try to find such a solution by picking appropriate initial and 
boundary values, but it is complicated by the large amount of initial and boundary 
parameters we have to choose. 

• For numerical purposes we have to introduce a second boundary at large m, which 
can introduce artificial boundary effects. 

• The initial values should represent a pre-classical function such that differences ap- 
pearing in the evolution equation are small. This requires a large lattice in the 
m-direction which increases the numerical effort to evolve and to solve the system of 
linear equations in order to find the initial values at the second n-slice. Furthermore, 
more m-values increase the likelihood of rounding errors adding up. 

• Close to classical singularities the wave packet can move rapidly on the lattice (corre- 
sponding to large extrinsic curvature). If this happens even between two consecutive 
n-slices, differences in the n-direction will become large which can lead to contribu- 
tions increasing exponentially toward the classical singularity. While imposing the 
consistency conditions exactly would eliminate those contributions, they will always 
be present in a numerical analysis and can even lead to wrong initial values due to 
rounding errors. 

In order to show an illustration of a possible pre-classical solution, Figs. Q and El have 
been obtained as follows: First, the evolution equation shows that for large n we have 
t n +i,m — tn-i,m, provided that the m-differences of t are small. Therefore, we can suppress 
possible oscillations in the n-direction, coming from consistent initial values at the second 
n-slice not close to the one on the first slice, by plotting only every second n-value. Secondly, 
in order to avoid a rapid movement of the wave packet in the m-direction, the lattice has 
been stretched in the n-direction by using the Hamiltonian Hk = ApA + 2kcp c . While the 
original value k = 1 leads to a rapid movement at small n which implies rounding errors 
due to exponentially increasing contributions and an unacceptable numerical solution, this 
turns out not to pose a problem for k > 4 (k = 4 in the figures) since the movement would 
happen between n = and n = 1 which is not visible. 

One could think that it is even impossible to find a pre-classical wave function reflecting 
the fact that De Witt's initial condition, which is close to the result of the dynamical initial 
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condition, is ill-posed in this model (i.e., its only solution would be zero). This, however, 
is not the case for two reasons: first, the notion of pre-classicality allows some tolerance 
since a discrete wave function can only be required to be close to a smooth one, but not 
exactly smooth; and secondly, the discreteness leads to a different behavior of solutions at 
small n such that a DeWitt-like initial condition is in fact well-posed (see [TT|). 
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